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Abstract 

We present a generalization of the coordinate Bethe ansatz that allows us 
to solve integrable open XXZ and ASEP models with non-diagonal boundary 
matrices, provided their parameters obey some relations. These relations 
extend the ones already known in the literature in the context of algebraic or 
functional Bethe ansatz. The eigenvectors are represented as sums over cosets 
of the BC n Weyl group. 
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1 Introduction 



In the framework of integrable systems, the study of open spin chains with integrable boundaries 
have been developed a long time ago [TJ [2]. There, it has been shown that the model is integrable 
provided the two matrices characterizing the boundaries obey some algebraic relations (the so-called 
reflection equation). However, although the integrability has been proven, the explicit resolution 
(eigenvalues and eigenvectors of the Hamiltonian) of the models is not known in its full generality 
In fact, for a long time, only the case of diagonal boundary matrices was solved, for different kind of 
models, e.g. open XXX [3] and its su(N) |4J or su(N\M) [51 [6] generalizations, open XXZ |7J and 
its generalizations [HI EJ [Hll EH], using different versions of the Bethe ansatz (analytical, algebraic 
or functional). However, the classifications of boundary matrices (obeying a reflection equation) 
[T2l [TBI [HI [15] clearly shows that non-diagonal solutions do exist, although explicit solutions for the 
eigenvalue problem were not known. The problem laid essentially in the construction of a reference 
state (a particular Hamiltonian eigenvector) that allows to initiate the procedure. Indeed, when the 
boundary matrices were not diagonal (or at least not simultaneously diagonalizable), the existence 
of this reference state was not ensured. 

Recently, different approaches have been developed to overcome these difficulties, such as gauge 
transformations that allow to go to a diagonal basis [16J, or fusion relations for TQ relations that 
do not need the existence of a reference state [TTJ [18]. In all these cases, the boundary matrices 
can be non-diagonal, but the parameters entering their definition need to satisfy some constraints. 
Note also the original approach [19] that uses another presentation of the reflection algebra, called 
g-Dolan-Grady relation, as well as the 'generalized' functional ansatz developed in [20J. Both avoid 
the use of constraints. 

The Asymmetric Simple Exclusion Process (ASEP) is an out-of-equilibrium statistical physics 
representation of the Temperley-Lieb algebra [21 J, on which the XXZ Hamiltonian is based. Many 
exact probabilistic results [22j [231 1211 [25], that are not necessarily based on the integrability of the 
model, have been obtained in the past fifteen years and shed some new light on integrable results. 
The ASEP notations, as described in [26J , are presented in this paper and a specific section 15.11 is 
dedicated to the comparison of our results with old results on the ASEP, such as the matrix ansatz 
[221 (23 [2H1 129]. 

In this paper, we present a construction based on the coordinate Bethe ansatz [30] for open XXZ 
and ASEP models, allowing the use of non-diagonal boundary matrices. As for the other approaches, 
the boundary matrices need to obey some constraints, but the ones we find are more general than 
those already known. 

We present now the structure of the paper: section [2] defines the models and introduces the 
notations needed for the next sections; section [3] describes the different sets of reference states and 
the sets of exceptional points they lead to. Section @] gives the structure of the coordinate Bethe 
ansatz as a combination of the reference states and explains the role of the BC n Weyl group and 
gives in details the Bethe equations that correspond to the two sets of specific points. Finally, we 
conclude on two open problems: section 15.11 presents a short discussion of the relation between the 
specific points and other previous results about the exclusion process such as the matrix ansatz; a 
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discussion on the completeness of the spectrum and of the eigenvectors is tackled in section 15.21 



2 XXZ and ASEP models with non-diagonal boundaries 



The Markov transition matrix for the open ASEP model is given by 

L-1 

W = K 1 + K L + y %2w jJ+1 , 
where the indices indicate the spaces in which the following matrices act non trivially 
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It is well-established [3T] 128] 126] that this ASEP model is related by a similarity transformation 
to the following integrable open XXZ model 

L-1 



where a are the usual Pauli matrices and 
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The explicit form of the transformation can be found in e.g. [26]. We reproduce it here for complete- 
ness: 
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where ^ is an arbitrary (gauge) parameter that disappears in all the following computations. 

Throughout the paper we will stick to this ASEP notation, keeping in mind that by this similarity 
transformation we treat also the XXZ model with non-diagonal boundaries. However, in the ASEP 
model, all the parameters must be positive (transition rates) while there is not such constraint in 
the XXZ spin chain. To be as general as possible, we will not assume this constraint in the paper. 

The vectors of the canonical basis used for the previous definitions are indexed by the spin value 
t or 4- for the XXZ-spin chain and corresponds to the number of particle Tj G {0, 1} of each site in 
the ASEP (see figured]). 
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Figure 1: Asymmetric exclusion process for a system of size L with two boundaries and mapping to 
the XXZ spin chain. 

3 Basis vectors for the coordinate Bethe ansatz 

For a chain of L sites with periodic or diagonal boundary conditions, it is easy to find one eigenstate, 
called pseudo- vacuum. It is usually chosen as L spins up (resp. L empty sites for the ASEP). Then, 
excited states are constructed by adding some excitations that flip a given number n of spins (resp. 
add n particles in the ASEP). The conventional Bethe methods (coordinate, algebraic or analytical) 

allow us to compute the linear combinations between the ( ^ j excited states that diagonalize the 

Hamiltonian. 

For an open chain with non-diagonal boundary (which is the case we want to deal with, see 
equations (12. lj) and (12. 3ft ). the first problem consists in finding the pseudo- vacuum. Different ap- 
proaches have been elaborated to overcome this problem as explained in the introduction. However, 
the coordinate Bethe ansatz, which is the historical method introduced by H. Bethe [HD], has never 
been successfully applied. To treat this problem, we must generalize the ansatz. A first step in this 
direction have been done in [21] . The general strategy we adopt here can be summarized as follows: 



• we do not choose anymore all spins up for the vacuum and spins down for the excitations. We 
take general vectors (see their explicit construction below), 

• Hamiltonian eigenfunctions will be constructed as linear combinations between states with n 
excitations together with m (< n) excitations, 

• the explicit forms of the excitations and of the vacuum depend on the total number of excitations 
in the state we consider. 

More physically, one allows that excitations may be destroyed or created by the boundaries with 
the restriction that one cannot create more than n excitations. Without any loss of generality, we 
can always choose that one of the boundaries preserves the number of excitations. In this paper, we 
choose the right-hand-side one. 

Finally, we found two different ways to fulfill these conditions. In both cases, some constraints 
between the parameters of the model appear (as it is already the case in the approaches [16j [26] 
to non-diagonal boundaries). We will show that we can solve the problem (within the framework 
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of coordinate Bethe ansatz) for all the sets of constraints that have already been produced in the 
literature as well as for some new sets. 

The following subsection introduces some notations that are standard for the ASEP and that will 
be used throughout the paper. Then, in the next two subsections 13.21 and I3.3[ we will present the 
two different sets of vectors with the associated constraints. Finally, in section HJ we present their 
particular linear combinations diagonalizing the Hamiltonian. 



3.1 Boundary operators and a duality 



Two particular sets of 2-dimensional vectors are relevant in the study of the boudary dynamics. The 
first set corresponds to vectors that diagonalize the boundary operators. We choose the following 
arbitrary normalizations, that will become clearer later: 
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(3.1) 



There are two solutions for the first equation (resp. the second) given by c\ = c±(a, r y) and Ai = 
A±(a, 7) (resp. cl = c±(j3,6) and A^ = A±(/3, 8)) where c±(u,v) and X±(u,v) are the roots of the 
two functions: 

P U>V (X) = uX + {u — v) — v /X i.e. c + (w, v ) = v/u and c_(w, v) = — 1 , (3.2) 
Qu,v(X) = X 2 + X{u + v) i.e. X+(u, v) = and A_(u, v ) = — u — v . (3.3) 

These notations allow for a convenient way of parametrizing both boundaries at the same time. 
The second set of relevant vectors satisfies the diagonal relations: 
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(3.4) 



As previously, the coefficients c* and X{ (resp. c* L and X* L ) take the two possible values c?j_(a,7) and 
A?j_(a,7) (resp. c±((3,5) and X±((3,5)), which are the respective zeroes of 



P n %( X) = uX + (u - v + q - p) - v/X = P UiV {X) + (q-p) 
Ql,vi X ) = X 2 + X(u + v +p + q) + {qp + pu + qv) . 

The explicit values of c* ± (u, v) and X* ± (u, v) are given by: 

p — q + v — m± a/ (p — q + v — u) 2 + Auv 



(3.5) 
(3.6) 



c*±(u,v) 
X*±(u,v) 



2u 



—p — q — v — u ± a/(p — q + v — u) 2 + Auv 



(3.7) 
(3.8) 



1 These functions are sometimes called k±(u,v). We change this notation to be consistent throughout the paper. 
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The definitions ( 13. 4 j) suggest the definition of the new operators: 
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where / is the two by two identity matrix. One verifies that the matrix W given in (12. ip has the 
following second representation: 

L-l 



w 



(3.10) 
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This representation will be essential in the definitions of the new exceptional points. 
3.2 First choice 

We first consider the matrix W as written in eq. ( 12. ip . Let us define the family of vectors: 



\u(u))i 



\V)i 



u{p/q) 

q-p 




i-l > 



(3.11) 
(3.12) 



where u is still an arbitrary parameter. From these elementary vectors, we build the tensor product 
over the sites i to j and write it as 



\fl{u))l = \u(u))i\u(u)) i+ i . . . \u(u))j . 



(3.13) 



We now fix an integer n and introduce the state with n — m excitations at the ordered positions 
1 < x m+i < ■ ■ ■ < x n < L defined as the (C 2 )® L - vector: 

\x m+1 ,...,x n ) =(./-) x (3.14) 



x |fi(w m+1 ))i m+1 1 \u)(v m+l )) Xm+1 \Q(u m+2 ))l2t 2 1+ i \u{y m+2 )) Xm+2 ■ ■ ■ \uj(v n )) Xn \VL{u n+1 )) L Xn+1 . 
The overall factor (w -) is introduced only in order to normalize the Bethe roots. The coefficients 



u m and related through the recursion relation: 

Q 



u m +i — —u m , l>m+l — — v m , (3.15) 

p p 

and the initial coefficients u\ and v\ are still arbitrary. These vectors correspond to states where m 
excitations have left the system (through the left boundary). To clarify the notation, the state with 
no excitation corresponds to m = n and is given by 

|0) = \n(un+i))t, (3-16) 
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while the state with one excitation (m = n — 1) reads 



\xn) = I^K))?'"- 1 Hv n )) Xn Mu n+1 )) L Xn+1 . (3.17) 

The above states are product states and are linearly independent. We denote by B n the set of 
these ( J independent vectors. We also need the vectors |x m +i, . . . ,x a , . . . , x n ) deduced from 

fc=0 ^ ' 

|x m+ i, . . . ,x n ) by replacing the vector \uj(u)) in position x a by |V)a; Q . 

These vectors have been chosen such that we get for the Hamiltonian bulk part: 



Xn) 



fO if/3,/3 + 1^ 

Xm+1 j • • • i Xn i 

y/pq\~,Xj - 1, ••) - q\..,Xj, ..) + \..,Tj, ..) H0 + l = x j ,p^ x S -x ; 

y/pq\..,Xj + 1, ..) -p\..,xj, ..) - \..,Xj, ..) if = xj , /3 + 1 ^ x j+1 ; 

v. | ••) Xji •Ej+i) ■■) I • •) Xj+i > ••) if /3 , /3 -(- 1 Xj+i. 

(3.18) 

Remark the additional states appearing on the right hand side involving vector \V). Because of the 
alternating signs, they cancel each other in the Hamiltonian bulk part, but the first and last terms 
(boundary terms). It is interesting to remark that the same "telescopic" trick is also used in the 
proof of the matrix ansatz in section [57X1 and in the definitions ( 13. 9ft . 

The first and last terms — |V)i and \V)l of the telescopic sum must be absorbed by the boundary 
operators and this justifies the study of the operator (13. 4ft . 

We now require that the left boundary diagonalizes \u(u\))i and that the right boundary di- 
agonalizes \u(u n+ i)) l, hence excitations are created neither on the right, nor on the left, when n 
particles are already in the bulk. Thus, the action of the boundary operators must be 



K\ \Xm+l i • • • i Xn) 



A-i l^m+l j • • • j Xn) 4" C\ 1 1, X m J r \ 1 . . . , X n ) if X m +i > 1 

A-l 1 1, Xm+2 •••) X n ) + D\ |^ m +2) • • • j X n ) + |1, S m +2) •■) X n ) if X m +1 = 1 



(3.19) 

K L \x m+1 ,...,x n ) = h l * m+1 -- Xn) _ l Xn< _ L > (3.20) 

I Aipm+i, • • • , X n -i, L) — |X m +i, . . . , X n -i, L) II X n — L , 

with the constraint 

Cf ] = . (3.21) 

We introduce also the more compact form Ai = A^ since this value will appear in the energy. 

The constraints (I3.2ip and (13.201) lead to the following trivial identifications with the parameters 
introduced in section 13.11 

ui = e s /c e (a, 7) , Ai = A e (a, 7) , (3.22) 

{p/q) L ' l u n+ i = c e ,(P, 6) , A L = \ e ,(/3, 6) , (3.23) 

(p/qf-'vn = c*Af3, 6), A L = \*Af3, 5)+q, (3.24) 



6 



where e, e', e" G {+, — } are arbitrary. The compatibility relation of the previous equalities with the 
recursion relations ( 13.150 lead to four possible solutions summarized in Table [TJ They correspond to 
the different choices of the signs e, e' in the relation: 



c € (a,j)c e ,(/3,5) 



L-l-n 



(3.25) 



The third sign e" does not appear in this condition, nor in the energy, nor in the Bethe equations: it 
just corresponds to the full reflection of the excitations on the right boundarjJE 



Ai 






a. 


7) 


c e ,((3,5) 


Constraints 








c+ 


> 


7) 




(e\ L 1 

7<5 \q J 


= 1 


—a — 7 


-(3-5 


C- 


[a 


7) 






1 


—a — 7 





C- 


'« 


7) 


c + ((3,5) 


-f^)™ 


= 1 





-(3-5 


c+ 


[a 


7) 




/ \ L— 1— n 

(?) 


= 1 



Table 1: Different possible values for the parameters and the constraints imposed by ( 13.250 . 



We then obtain the following values of the coefficients: 

,(m) _ u 2 



(m-l) 



Vi - U 2 



V! - U 2 



-P , jJ (e s /u m +i) 



Ai 



(m) _ (g/p) m (Ai + a)vi - (p/g) m (Ax + j)u 2 + 7^2 - avi 



Vi - u 2 



Ai' 



(m-i) _ (p/?) m (Ai + 7)1*2 - (q/p) m v 1 (A 1 + a) - 7^1 + (a + g - p)w 2 



(3.26) 
(3.27) 

(3.28) 

(3.29) 



We also remind that for ASEP probabilistic models, all parameters have to be positive, so that 
only the two first lines of table [1] have to be considered in this case. These two constraints can be 
recasted into a single one 



n / \ L—l—n / \ L—l+n 



0. 



J) J / \ \p y 
Then, using the correspondence ( 12.80 . it takes the form for XXZ model for even L: 

cos(6 l i — 6 2 ) = cos(o; + + w_ — mrj) , 



(3.30) 



(3.31) 



that is just the constraint given in [TTJ [T6], [26] for some integer m. Doing the same with the last two 
constraints, we get 

cos(#i — 62) = cos(cj + — uj- — nrj) . (3.32) 



2 Of course, we could have also chosen a full reflection on the left boundary instead. 
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3.3 Second choice 

Using the same technique as in the previous subsection, we introduce a basis of states that is suitable 
to the diagonalization of the alternative representation (I3.10p of the matrix W. The vector \V)i has 
the same definition as before but we introduce the family of vectors 

\ U *{u))i = Q • (3.33) 

We define similarly the product state as in (13. 13[) by replacing \co(u)) by \cj*(u)) and we 

introduce the new state with the new excitations in ordered positions 1 < y m +i < ■ ■ ■ < Vn < n\ 

Iw, ■ ■ .,y n y = ( A /|V m+1 1)+ " My " ^ l^*« l+ i))r +1 "V*K + i)), m+1 • • • KK))yJ^(< + i)>L 



(3.34) 



9, 

where the new coefficients u* and v* satisfy the recursion relation 

Kn+l = > V *m+1 = (P/V) v m ■ ( 3 - 35 ) 

These vectors have been chosen such that we get for the Hamiltonian bulk part 



* I \ * 



fO if/3,/3 + 1^ 

- !••)* - P\ -Vr)*- I ••!£••)* if/9 + l = %- ; P^Vj-i] 
pq\..yj + l..)* -q\.., yj,. .)* + {.., yj, ..)* if = yj ; (3 + l^y j+1 ; 
P, ..)* - \..,J+T, ..)* if P = ; /3 + 1 = . 

(3.36) 



As for the previous choice, there are additional states |V)i and \V)l that survive the telescopic 



sum J2i=i w ii+i an d that must be absorbed by the boundary operators. We impose the dynamics: 



A*.H| r \*_L/7* ,(m) M T , \* if T* 1 ^> 1 

Km+li • • •/ i ^ l I -L? x m+l) • • •/ 11 x m+l x 

A^ ^ |1, x m+ 2 ■■■)* + ^|x m+ 2, ...)* — |1, X m+ 2, ... )* if X m+ \ = 1 



jvr*. \* _ I A2|x m+ i, . . . ,x„)* ifx n <L, 

with the additional closure constraint 

C*' (0) = . (3.39) 

We introduce again the more compact form A^ = A{ since it is this value that will appear in 

the energy. From the constraints ( I3.38P and (13.391) and the notations introduced in section 13. H we 
identify directly the value of most parameters: 

uj = e7c:(a, 7 ), A* = A:(a, 7 ), (3.40) 

< +1 = c* t/ ((3,5), A* L = \*A(3,5)+ P + q (3.41) 

< = c e »CM, A* L = p + \ e// (P,8). (3.42) 
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As in the previous subsection, these values are compatible with the recursion relation f l3.35[) if and 
only if the parameters of the model satisfy the relation 



c:(a, 7)408, *) = 




(3.43) 



for some signs e and e' and some integer n, which then fixes the number of excitations. Similar 
exceptional points already appeared in the literature in the context of the matrix ansatz described 
in section 15.11 Table [2] summarizes the different possibilities. 
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c;(a, 7 )c;(/3,5)=e* (j) 
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-a-y-P-S-y/ (p-g+7-a) 2 +47Q- v /(p-g+5-/3) 2 +4 / 3<5 






clK 7 )cl(/3,5)=e 5 (V) 
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-a-j-p-S+y/ (p-g+7-a) 2 +47a- v /(p-,j+(5-/3) 2 +4/3<5 
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c;(a, 7 )c!L(/3,5)=e 5 (|) 
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-a-i-p-S-y/ (p-g+7-a) 2 +47Q+ v / (p- y +(5-/3) 2 +4/3<5 
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cl(«, 7 )c;(/3,5) = e s (j) 
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2 



Table 2: Different possible values for the parameters and the constraints imposed by Q3.25]) . 



The values of all the coefficients are then given by: 

c *(m) = P; 7 (eV< + i) ( 4) 

1 2 v\ - u* 2 y ' 

D .,(m-1) = _ v , P a „{eyv*J ^ (3 45) 

A *,(m) = (p/qrvKAl + a + q)- (q/p) m u*(Al + 7 + p) + (7 + ~ (<* + gK ,3 4g) 

^*,(m-i) = (g/p) m ^ 2 (A^ + 7 + p) - (p/g) m (A* + a + gK + (g + P>*2 - (l + PM ,3 4?) 

In the language of the XXZ spin chain, we should introduce the additional angles Q± defined by 

2 cos Q± = cos S± — cos u± , (3.48) 

and the condition (I3.43[) takes the following simple form for even L: 

cos(#i - 6 2 ) = cos{tt- ± 0+ - rjm) , (3.49) 

for some even integer m between — L and L. Up to our knowledge, this condition is new. The 
associated Bethe equations are presented in section 14.21 
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4 Coordinate Bethe Ansatz 



4.1 The first set of specific points 

We are now in position to propose an ansatz for the eigenfunction 

$ « = E E £ 4 w) ^ Cm) iwi,. -,*„>, (4.i) 

m=0 x m+ i< - <x„ geGm 

where G m is a full set of representatives of the following coset BC n /BC m (see appendix [SJ, the 
vectors \x m+ i, . . . , x n ) are either given by ( 13.141) or by ( 13.341) and we introduce the notation for 
the following truncated vector 

k^ = (k m+1 ,...,k n ). (4.2) 

For this definition to be consistent, the coefficients A g m ^ do not have to depend on the choice of the 
representative i.e. 

A^ = A^ for any h G BC m . (4.3) 

The coefficients A g m ^ are complex numbers to be determined such that $„ is an eigenfunction of 
H i.e. such that the following equation holds 

H® n = E§ n . (4.4) 

Due to the results of section [31 HB n C B n . Then, we project equation (14.41) on the different indepen- 
dent vectors of B n to get constraints on the coefficients A g m \ This projection depends on the two 
choices done for the vectors |x m+ i, . . . ,x n ). However, we chose the coefficients in such a way that 
one can be deduced from the other easily by adding a star to most quantities. Therefore, we write 
only the projections for the choice of section f3T2l and only the ones leading to independent relations 
(one can check that the remaining ones do not lead to new relations). 

On \xi, . . . , x n ) for (aii, . . . , x n ) generic (i.e. 1 < x%, x n < L and 1 + Xj < Xj + i) 

As in the usual coordinate Bethe ansatz [30], this projection provides the energy: 

E = Ai+A L +^ A(e ifc where \{x) = ^pq(x+^-) -p-q = ^ (x - (x - . (4.5) 

Let us remark that, up to the boundary terms Ai and A^, the energy takes the same form as in the 
periodic case. 

On \xi, . . . , x n ) with Xj+i = 1 + Xj (and x\, . . . , Xj-i, Xj + 2, ■ ■ ■ , x n generic) 

This projection is also a usual one and provides the scattering matrix between excitations. It is 
given by a relation between Ag and A g °J j where <jj is the permutation of j and j + 1 (see appendix 
|A|) . Namely, we get 

A<®. = S (e ik °\e ik ^+v) Af , (4.6) 
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with 

S(u,v) = w here a(x,y) = — % — ( ( ] y - xy - 1 ] . (4.7) 

a(v,u) xy - 1 V VV P V 9/ / 

The normalisation chosen for the function a(x, y) is for further simplifications. As expected, this 

relation is similar to the periodic case since the boundaries are not involved in this process. 

On |1, x m+1 . . . , x n ) (x m+1 , ...,x n generic and m > 1) 

This relation is a new one and we must take into account that the left boundary can create one 
particle at the site 1. We finally get, for any g G G m , 

m 

(l^ -Ai-P + J2(p + < 1- VP~q{e lkaj + e"**"-))) Yl A t~ 1)elkahm 

3=1 h£Hm 

+VM Yl A^- 1] e 2ik ^ + C[ m) A g m ^ = , (4.8) 

heH m 

where H m = BC m / BC m -\. To obtain this relation, we have used the following property 

J2 A g e ik ^e ik ^ )x(m) = J2 e ik v m)x{m) A gh e ik ^ . (4.9) 

g&Gm-i g&Gm h£H m 

Let us stress that equation (14. 8 p is invariant by the choice of representative of G m . 

On \x m+1 . . . , x n ) (x m+ i, ...,x n generic and m > 1) 

This projection provides a second relation between the coefficient from the level m — 1 and m 
since we must take into account that the left boundary can destroy a particle present on the site 1. 
We obtain the following constraint, for any g G G m , 

m 

D (m-i) Aj^'e^w + (A[ m ) _ Ax + Y,(P + V - VPV( eiksJ + e-^)))4 m) = . (4.10) 

heHm j=l 

From (14. 8p and (I4.10p . we may express all the A g m ^ (m > 1) in terms of A g thanks to the following 
recursive relations defined for m > 1 



A (rn) = T M( e ifc sl) _ _ . ^ k 9^ A {rn-l) , (4.11) 



with the following definitions: 



\xi, . . . , x m ) — - - m _ 1 - - - — - , (4.12) 

Pi[x m )Vi[x m ) l[. =1 a(x m ,x i )a(x i , l/x m ) 

V X {Z) = X(Z) + (Ax + 7 ) (l - ^) + (Ai + a)(l - , (4.13) 

p 1 (Z) = Z + r, (4.14) 
1 pu 2 -qv 1 1 pu m+1 -qv m 

(4.15) 
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where r describes the difference between \u>(vk)) and the vacuum states \u(uk)) and \u(uk+x))- It is 
independent from k. 

The proof that (14.111) is a solution of both equations (14 .8p and (14.101) is postponed to appendix [B] 
and relies on a residue computation. The integrability of the model plays a role at this place, since 
there are a priori too many constraints but not all of them are independent. 

Let us emphasize that this recurrence relation (14.111) is the main result of this article. Indeed, 
firstly, it proves that the ansatz (14.11) and, in particular, the choice of cosets used to define it, is the 
appropriate choice. Secondly, it gives a very simple way to construct eigenfunctions, and we believe 
that it may be used for further computations. 

A consequence of ( 14. lip (for m = 1) and Ag\ = A g , is 



A (o) = I^± A (o) . 



"9ri / e -ik B iy ~9 ■ (4.16) 

This relation, together with ( 14.61) . allow us to express A g ^ for any g G BC n in terms of A^ (where 
the subscript 1 stands for the unit of BC n group). Finally, using recursively (14. lip , we can express 
all the coefficients A g m ^ in terms of only A^. This last coefficient is usually chosen such that the 
eigenfunction $ n be normed. 



On |a?i . . . , :c n _i, L) (x 1 . . . , x n _ x generic) 



This last constraint consists in the quantization of the excitations moments since the system is 
in a finite volume. In the context of the coordinate Bethe ansatz, this quantization leads to the 
so-called Bethe equations, explicitly given by, for 1 < j < n, 

FT S(e* e^Sie-**^**) = e 2lLk > ' [ ik ' L ) ' , (4.17) 

V L (x) = X(x) + (A L +P)(l- + (Al + . (4.18) 

V x y p/ V x y qJ 

We remind that the parameters have to obey one of the relations given in table [U 



4.2 The new second set of specific points 

When condition ( 13.431) is satisfied, the same construction as in the previous section works, up to easy 
modifications in the equations. We introduce a new function V±{Z) and a new parameter r*, which 
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replaces r, defined by: 



V*(Z) = \(Z) + (Al + a + q) II- 
Vl(Z) = X(Z) + (A* L + a + q ) (l - 

1* * 



1 /J) 

1 /p 



+ (A* + 7 +p) ( 1 — - 
+ (A2 + 7 + p) f 1- 





(4.19) 


V p/ 




(4.20) 


VP/ ' 




(4.21) 



y/pq u*-vl 

The same proof as the one presented in appendix |A] is valid and we obtain the Bethe equations for 
3 e {l,...,n}: 

e H \ . , u = TT^(e ifci ,e^)5(e-^,e ifci ) , (4.22) 

with the parameters obeying now one of the relations described in table [2j 



5 Conclusions 

We discuss in this section two open problems, which need further investigation. 



5.1 Connection with the matrix ansatz 

For s — 0, the matrix W is stochastic (Markov transition matrix of the exclusion process) and hence 
has a known ground state eigenvalue E = with a simple left eigenvector, whose components are 
all equal to 1 (conservation of total probability). Nevertheless, the corresponding right eigenvector 
is non-trivial since the matrix is not Hermitian and describes the stationary properties of the asym- 
metric exclusion process. The structure of this specific eigenvector was elucidated first in [22J and 
then studied algebraically in more details in [2"gj |2"T]. 

As already explained, the state space is 2 L dimensional and the canonical basis can be indexed 
by the values of the occupation Tj G {0, 1} (resp. spin Sj G {—1, 1} in the XXZ language) on each 
site i. The matrix ansatz states that the ground state of the ASEP with s = has components given 
by: ^ 

(r 1 r 2 ...r L \^) = ((V 1 \ J] ( Tl D + (1 - n)E) \V 2 }} , (5.1) 

l<i<n 

where the arrow means that the product have to be build from left to right when the index i increases. 
One has for example (00 . . . 00|$) = ((V^E 1 ^)). The non-commuting matrices D and E act on an 
abstract auxiliary vector space V. The vector |Vi)) is in this space V, whereas the vector ((V 2 | is in 
its dual. $ is an eigenvector of W for E = when s — if the two matrices D et E and the two 
boundary vectors satisfy the commutation rules [22J: 

pDE - qED = D + E, (5.2a) 
((V 1 \( 1 D-aE) = -((V 1 \, (5.2b) 
((3D-5E)\V 2 )) = \V 2 )). (5.2c) 



13 



with the condition s = 0. These three relations allow to determine recursively all the components 
of the eigenvector $ and do not need an explicit representation of the algebra. In particular, the 
matrix ansatz gives an easy access to classical correlation functions with standard transfer matrix 
techniques. 

The connection with the present Bethe ansatz approach comes from the algebraic study of the 
algebra generated by D and E, as performed in [281 [27], from computations that arise in the combina- 
torics of some so-called staircase tableaux [29] , as well as from the physical interpretation [23l [25l [24] 
of the matrix ansatz in terms of so-called "shock" product state as introduced in section [3j 

The first remark is that the matrix ansatz described here fails for some values of the parameters 
p,q,a, {3,j,8. Failure happens when the recursion relation on the size L induced by (15.21) leads to 
((V2IV1)) = and thus more generally to a null vector $, as explained in |28j. This failure condition 
for a system of size L is the existence of an integer n & {0, 1, . . . ,L — 1} such that 



This relation is precisely one of the two cases of ( 13 .30 j) when s = 0, for which the present Bethe 
ansatz approach works. 

The second remark is that the operators D and E do not have generically finite dimensional 
representations. Detailed studies shows that such finite-dimensional representations may exist at 
some specific points [271 EE]- It appears that an ra-dimensional representation exists if one among 
the four cases of the condition (I3.43|) is satisfied for some integer n and plus/minus signs, i.e. 



In [27], the finite representation is used to study how the ground state evolves through the phase 
transition that occurs in the open ASEP along the manifold (15.41) : the present paper gives also the 
excited states on the same manifold and the Bethe ansatz equations we obtain in this case now allow 
one to study how the gap behaves near this phase transition. Once again, atypical results appear in 
both the Bethe ansatz approach and the matrix ansatz approach for the same specific parameters. 

A third remark relies on the observation that the present construction does not allow to express 
the ground state described by the matrix ansatz in terms of a coordinate Bethe ansatz, precisely 
because one fails when the second works. This mismatch leads us to think that the matrix ansatz 
state may be used as a new reference state to build the missing eigenvectors if one could manage to 
add excitations on it. However, up to our knowledge, no deep understanding of the relation between 
these two approaches exists and has been exploited. 

A last remark deals with the combinatorics work [29]. Besides the fact that quantities such 
as (I3.30p and (I3.43j) appear in numerators or denominators in their computations and thus either 
simplifies or invalidates the matrix ansatz, their rewriting of the matrix ansatz in terms of staircase 
tableaux sheds a new light on the precise structure of this ansatz, and puts it in a form that is more 
suitable to comparison with Bethe ansatz because of the relation of these tableaux with permutation 
tableaux. In particular, the number of so-called staircase tableaux is n\ 4 n , whereas the cardinal of 




(5.3) 




(5.4) 
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BC n is precisely n\ 2 n . Understanding both the similarities and the mismatch may lead to a better 
understanding of the Bethe ansatz. 

5.2 Completeness of the spectrum 

As usual for Bethe Ansatz methods, the delicate point to check is the completeness of the spectrum. 
Indeed, in many cases, completeness is not proved but expected to be true: numerics for small size 
systems [32], enumeration of the number of roots in the thermodynamic limit |33j. 

In the present case, numerical checks of the completeness for small sizes have been considered in 
|34j : the full spectrum is shown to be described by two sets of Bethe equations. Our approach for 
the right eigenvectors gives only one set of Bethe equations for each specific point and thus only one 
part of the spectrum and one part of the right eigenvectors. The same construction can be done for 
left eigenvectors (or right eigenvectors of the adjoint). In this case, one obtains the other part of the 
spectrum and one part of the left eigenvectors, as discussed in [21] . 

Let us stress that the operators we are considering may not be hermitian and thus left and right 
eigenvectors may not be simply related. Although the full spectrum is known, it would be interesting 
to build the full set of right and left eigenvectors. 
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A Weyl group BC n and coset 

The Weyl group BC n is generated by the set {r\, o"i, . . . , <J n -i} with the following constraints: 

cr| = 1 = r\ , a 1 r 1 cr 1 r 1 = r 1 aina 1 , o/r, . yt j j = a j+1 aja j+1 . (A.l) 

It acts on a vector k = (ki, . . . , k n ) of C n : 

k ri = (— ki, &2, • • • > k n ) , k a . = [k\, . . . , fcj+i, kj, . . . , k n ) . (A.2) 

The subgroup generated by {<7i, . . . , 0Vi-i} is just the symmetric group. We now consider its sub- 
groups generated by {r l5 <7i, . . . , <r m _i}, m < n, which we identify with BC m . 

For g G BC n , we then define the class gBC m = {gh; h e BC m }, called a left coset. It is known 
that the set of all classes gBC„, which is called BC n /BC m , forms a partition of BC n : we can thus 
define G m as a full set of representatives of BC n / BC m , such that one has the unique decomposition 
BC n = © 9gGm gBC m . We set, by convention, Go ~ BC n and G n = {1}. 

The action of an element g of G m on a vector k^ = (k m+ i, . . . , k n ) of C n ~ m is given by 
k^ = (A; 9 ( m +i), kg ( n )). One checks that this action does not depend on the choice of the 
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representative g, such that the action of BC n / BC m is well-defined on C n ~ m without further spec- 
ifications. This definition is useful because the set {ky^\g e G m } contains one and only once 
the vector (e^k^, . . . , e in _ m k in ^ m ) for any choice ej — ±, 1 < ij < n and ij ^ ik- For example, 

{k<?~% e G n -i} = {(4)7 (-*„), (K-i), (-K-i), • • • , (-h)}. 

Finally, we introduce H m which is a full set of representatives of the coset BC m j BC m _\ used in 
section H] and in the next appendix to simplify computations. 



B Proof of solution (14.111) 

We want to prove that ( HTTTj) implies (J48]) and fTCTOj) . 



We start with ( I4.10p . We start by remarking that a consequence of (14. lip is 



1 < a < m — 1 , 



A {rn) = A (rn) x J (e^ e^- ) ( , fc jfc9m+l) 

'/fT,i i/ i Tl m ) (e s 1 ... e 9 m ) 



S(e 



« > m + 1 . 



(B.l) 



Then, using again ( 14. lip to express now Ag +1 ^ in terms of Ag U> and using ( IB.ip to express A { ^ > in 
terms of A^ , relation ( I4.10p becomes the functional relation 



,(m) 



,(m) 



m+1 

E 



m+1 



m+1 



X. 



Vi(xj)pi(xj) Y[ xe)a(x e , — ) + — Vi(— )pi(— ) J~J a(— , x^.)a(^, xj) - A(x 3 ) 



^=i 



ry* rr> ry* 

^3 3 ■ x '3 



OC i OC i 

3 e=1 3 



Ai - Ai 



(m+1) 



(B.2) 



where Xj stands for exp(ik g j) and the functions are defined in (I4.7p - (l4.13p . 
To prove this last relation (1B.2|) . let us introduce the following function 



F {m \Z) 



Vi(Z) Pl (Z) 



\{Z) (2Z-^ q + Jj 



m+1 ^ 

I [ a(Z,x e )a(x £ , —) 



(B.3) 



i=i 
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such that 

Res{F^ m \Z)) 



x. 



Z=Xi 



m+l j 
f — 1 



=1 



m+l 



Res(F^(Z)) = 1^(1)^(1) ]la(-,x e )a(x e 



Z=l/x 



Res(F (m) (Z)) 
Res(F (m) (Z)) 
Res(F< m >(Z)) 

Res(F (m) (Z)) 



v 7 ^ ( - 



m+l 



Ai + a ( fp 

— h r 



p — g \ y g 



p\ Ai + 7 / /g 
g/ p-q \\ P 



Z=oo 



2(p -_ (p _ ?+27 _ 2Q) Q(^ + ^ 1+rl , 



_g AM _^(i(^ + yD +r )_i (2Al+Q+7) 



(B.4) 
(B.5) 

(B.6) 
(B.7) 
(B.8) 

(B.9) 



Then, (1B.2j) is equivalent to ^ resi(iMe F^ m ^ (x) = which finishes the proof. 
For (14. 8p . we use the same procedure and (1B.2P to get 



m+l 

i>q Yl 

3=1 
m+l 



m+l 



m+l 



' 1111 X 

a;JVi(a; i )pi(a; j ) a(xj, x^)a(x£, — ) + — Vi(— - )pi(— ) J\ a(—,x e )a(xi, 



m+l 



J 3 3 J £ = i 3 

(m+l)\ ^-y(m+l) n (m) 



( ^ X( Xj ) + A x - AS m) + p) ( Hxj) + A i - A f n+1) ) " C[ m+1) D[ 

3=1 3=1 



(B.10) 



The function to consider is now G (m) (x) = ^/pqx 

_p( m )( x ). its residues at the point x = xq with 



xq 0, oo are simple to compute: Res(G ( - m \x)) 



pqx Res(F^ (x)) 



X=Xg 



. Since is not a 



pole, it remains to compute the residue at infinity. It reads: 



X=Xg 



m+l 



Res(G {m) (x)) 

x/pq 



y2X(x j )X(x k ) + 2[Res(F^\Z)) 

j,k=l v 



p + q 



Z=oo 



m+l 
3=1 



2 V r H\/f + yf))( 2A ' + a + ^ + ^ i )- (P "'4 a " 7> +1 H il! - (B11) 



Then, as above, (IB.lOp is equivalent to J2 residue G^ m \x) = 0. 
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